Large Magnetoresistance Induced by Quantum Charge Fluctuations in Magnetic 

Double Dots 
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We study electron tunneling through two small ferromagnetic dots. Quantum charge fluctuations 
and interdot coupling make each Coulomb peak of the conductance at zero interdot coupling split 
across. The interdot tunnel coupling is determined by the relative orientation of magnetizations of 
the two dots, leading to different splitting energies of the Coulomb peaks in parallel and antiparallel 
magnetization alignments. As a result, a very large tunneling magnetoresistance occurs near the 
Coulomb peaks, and its sign may be either positive or negative. 
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The tunneling conductance of a tunnel junction made 
of two ferromagnetic (FM) electrodes separated by a non- 
magnetic (NM) insulating layer may increase by several 
tens percent, as the magnetizations of the FM electrodes 
change from antiparallel (AP) to parallel (P) alignment 
under an external magnetic field (Ij, iS, iSj i^] . Such a 
tunneling magnetoresistance (TMR) effect has important 
applications in magnetic sensors, random access memo- 
ries, and magnetic imaging. The TMR originates from 
the spin polarization of tunneling electrons in the FM 
electrodes. The TMR ratio can be defined in a symmet- 
ric form as 



TMR: 



Gp — Ga 
Gp -l- Ga 



(1) 



where Gp and Gap are the conductances for parallel and 
antiparallel magnetization alignments, respectively. Ac- 
cording to the simple model P, 0, TMR = PlPr for 
the magnetic tunnel junctions, where Pl(r) is the spin 
polarization factor in the left (right) FM electrode. For 
electrodes with spin polarization of the same sign, the 
TMR ratio is usually positive. There is great interest in 
inverse (neg ative) TMR |ili,[ll[il[l3, 

which is usually 

associated with certain hiding mechanisms that affect the 
spin-dependent transport significantly. The inversion of 
TMR has been attributed to the reversion of the spin 
polarization in either of two FM electrodes due to the 
bonding effect at the FM/NM interface 8, 9, 10], and to 
the resonant tunneling via localized states on impurities 
in the insulating barrier 0, 0| . 

With the advancement of technology, magnetic tunnel 
junctions are made increasingly small and controllable in 
recent years. Single electron transistors consisting of FM 
metals ^3 and junctions with small FM islands have 
been fabricated. Presently, it is of particular interest to 
study the TMR through small junctions and quantum 
dots, where the Coulomb interaction plays an important 
role lis] , ft has been suggested that the 

TMR between two FM leads through an FM dot in the 



Coulomb blockade regime may be enhanced by higher- 
order electron tunneling via virtual intermediate states in 
the dot (cotunneling) IllTi llfilj or oscillate in magnitude 
with external voltage [13, lla]- The TMR through an 
NM dot may be unch ang ed or in some cases decreased 
by Coulomb blockade EE [3. 

Nonmagnetic quantum dots have previously been ex- 
tensively studied, and many features of the tunneling 
conductance have been well known The conduc- 

tance of a quantum dot displays Coulomb peaks at cer- 
tain gate voltages, and is strongly suppressed elsewhere 
by the Coulomb blockade jH [^l I3. Electron 
transport through coupled double quantum dots has also 
been extensively explored. In a double-dot structure, the 
quantum charge fluctuations can be more directly probed 
than in a single dot. As has been studied both experimen- 
tally Hi] and theoretically [H HI HI], the interdot cou- 
pling and charge fluctuations make each Coulomb peak 
at zero interdot couplings split across. However, spin- 
dependent electron transport through coupled FM dou- 
ble quantum dots has not yet been investigated so far. 
More importantly, we find that, in the FM double quan- 
tum dots, a largest TMR effect may occur essentially 
because of the quantum charge fluctuations between the 
FM dots. 

In this paper, we consider a magnetic tunnel junction 
consited of two coupled FM quantum dots. The sys- 
tem is weakly coupled to two reservoirs of electrons. A 
quantum-circuit description for the charge fluctuations 
is developed. We show that, with changing the magne- 
tizations of the two FM dots from AP to P alignment, 
the interdot tunnel coupling increases, which in turn en- 
hances the Coulomb peak splitting. As a consequence, 
the Coulomb peaks shift to different gate voltages for 
AP and P alignments, resulting in a very large TMR. 
Moreover, depending on the gate voltage, the TMR ratio 
may be either positive or negative. 

The equivalent electrostatic circuit for the double-dot 
tunnel junction is shown in Fig. 1. For simplicity, the 
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FIG. 1; Equivalent electrostatic circuit of a magnetic tunnel 
junction consisting of two small ferromagnetic quantum dots, 
which are weakly coupled to two electron reservoirs. The 
arrow represents the magnetization direction of the FM dot. 



two FM dots are assumed identical. For the symmetric 
capacitors: Cl = = Co and Cgi = Cg2 = Cg, we can 
confine ourselves to the case of the same gate voltage 
Vgi = Vg2 = Vq. The qualitative results of this work 
can be reproduced for nonidentical dots, provided that 
the two gate voltages V^^i and Vg2 are allowed to change 
differently so that electrons tunneling through the two 
dots resonates simultaneously. Let iVi and A'2 denote 
the numbers of excess electrons in the two dots. The 
model Hamiltonian for the system can be written in the 
form 



H = K +U{Ni,N2) 



(2) 



where the expression for electron kinetic energy K will be 
derived later, and the electrostatic energy of the system 
is given by 

U{Ni,N2) = Ec{6^Ni + S^N2) + 2E'JNiSN2 . (3) 

Here SNi = Ni - Nc and 6N2 = N2 - with A^c = 
CoVa/e. Ec = e^C^/2F and E'^ = e^C^j^F with = 
Co + Cg + Co and F = C^ - C^. 

The interdot tunnel coupling of electrons is assumed 
to be much stronger than the tunnel coupling between 
the dots and the reservoirs. As a result, for given total 
electron number M = N\ ^ N2 , quantum fluctuations of 
M — Ni — N2 are important. In previous treatments of 
the charge fluctuations, very different assumptions were 
made on the interdot tunnel matrix elements, e.g., from 
the thin-shell model with one-to-one level coupling 
to the thick-shell model with single site-to-site connec- 
tion I23, [13 ■ Here we extend the approach developed by 
Girvin et al. to the charge fluctuations in a Coulomb- 
blockaded single junction coupled to its environment by 
long transmission lines. In the approach, the collective 
motion of the charges is described by inductances, and 
the charge numbers are treated as quantum canonical 
coordinates. In this essence, we introduce an inductance 
I between the dots to take the inter-dot charge fluctu- 
ations into account. The advantage of the approach is 
that it enables us to capture the essential physics of the 
charge fluctuations without detailed assumptions about 



the tunnel matrix elements. It will be shown that a sim- 
ple connection can be established between our theory and 
the previous theories 113, . 

Since inductance I essentially describes the charge 
transfer between the dots, it must depend on the mag- 
netization alignment of the FM dots because of the spin 
polarization of electrons in the FM dots, just as in an 
ordinary FM/NM/FM tunnel junction For the P 

(AP) alignment, the electron tunneling between the dots 
is relatively easy (difficult), and so the inductance / = 1^ 
{I = Iap) is relatively small (large). According to classical 
electromagnetism, as a current of iVie — ~N2e — Aie/2 
flows through inductance with 77 = P or AP, the en- 
ergy stored in the inductance is K — /^e^TM^/S. Its 
corresponding quantum form is K = 2P^/l^e^ with Pj^ 
the momentum operator conjugate to ^A. Projecting K 
into the subspace of integer electron numbers A'^i and 
N2 113, we obtain K = -(£>!, + - 2)/2lr^e'^, where 
operator — g^^iP^/ri increases M. by two. Subse- 
quently, K is transformed into the representation of A'^i 
and N2, yielding 



K^-Kr,{D\D2+DiD. 



(4) 



where k,, — d\^^-^ — e *^i(2)/?i increases iVi(2) 

by unity with Pi{2) the momentum conjugate to ^^1(2), 
and an irrelevant constant energy l/ljje^ has been ne- 
glected. Equation (4) is a reasonable result, indicating 
the movement of an electron from one dot to the other. 

The above model is essentially a quantum LC circuit, 
whose general eigenstates can be determined from the 
complicated Mathieu equation |3nl |. In this work, un- 
der the assumption of Ec ^ E'^ and Eq 3> Hn and 
so to a good approximation, we project the Hamilto- 
nian into the subspace spanned by the lowest-energy 
states for all possible Af. The eigenstates of A''! and N2, 
\ni,n2) = (d1)^i (£'1)^2 10,0), can be used as the base 
wave vectors, where |0, 0) stands for the state without 
excess electrons in the dots. For even Af = 2k, the low- 
energy state is the single state \k,k) with electrostatic 
energy Ej^ = U{k, k); and for odd Af = 2k + 1, the low- 
energy states are the two degenerate states |fc -I- 1, fc) and 
\k,k + 1) with electrostatic energy Ej^ — U{k + l,k). 
Since the kinetic energy given by Eq. ^ conserves the 
total electron number Af, the only nonzero matrix ele- 
ments of K in the low-energy subspace are those between 
the two degenerate states for odd Af = 2fc -t- 1, namely, 
{k + 1, k\K\k, k + 1) = {k,k + l\K\k + 1, fc> = -k^. It 
is thus easy to obtain the eigenstates of the system. For 
an even electron number (Af = 2k), the eigenstate is 
|7V) = \k,k) with eigenenergy = Ej^. For an odd 
electron number {Af = 2k + 1), there are two eigenstates, 
which will be distinguished from each other by symbols t 
and i, symmetric state \N t) = i\k + l,k) + \k, k + l))/V2 



with lower eigenenergy Sj^-^ — Ej^ — and antisymmet- 
ric state \Af I) = {\k + l,k) - |fc, /c + l))/\/2 with higher 
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eigenenergy f^-i = Ej^ + Kj^. 

The total excess electron number M in the ground state 
is determined by the minimum of system eigenenergy. 
Therefore, tunneling of an electron into or out of the dots 
usually causes an increase of the system energy and so 
the tunneling rate is suppressed, leading to the Coulomb 
blockade effect. At certain values of the gate voltage, 
however, the system eigenenergies are degenerate for dif- 
ferent A/", i.e., E2k = E(2k+i)\ or E2k+2 = -E(2fe+i)|, in 
which there appear Coulomb peaks of linear conductance. 
The positions of the Coulomb peaks are obtained as 

V.^V4.^(. + i4) , (5) 

where = {E'^ + k^)/Ec describes the correction to the 
peak position due to interdot electrostatic coupling E'^ 
and interdot tunnel coupling k^. According to Eq. lO, 
if the interdot coupling is absent, the conductance peaks 
occur periodically, whenever the gate voltage is equal to 
half-integer times e/Cc- The presence of the interdot 
coupling splits each peak into two, the distance between 
them given by b^ejCc. It is interesting to notice that 
our result recovers the previous result for the NM dots 
obtained in the weak-tunneling regime [l^, 113 , if we set 
= (?iln2/7r)GD/Cj; with the interdot conductance. 
It then follows that, for the FM dots, does depend 
on the relative orientation of their magnetizations, for 
tunneling conductance Go between two FM dots must 
be orientation dependent. 

According to Eq. ISJ, there is a position shift of the 
Coulomb peaks with changing the magnetization align- 
ment from AP to P. This effect originates from the quan- 
tum charge fluctuations and results in a very large TMR. 
To evaluate the conductance and TMR, we will confine 
ourselves to the vicinities of the Coulomb peaks, where 
real tunneling processes of the electrons dominate. In the 
strongly Coulomb-blockaded region, higher-order tunnel- 
ing processes of electrons via virtual intermediate states 
(cotunneling) may be important [isl Il6l| , which will not 
be considered here. Let us consider the first pair of res- 
onant peaks for positive gate voltage Vgi whose posi- 
tions are determined from Eq. (0) by setting fc = 0. The 
physics around other peaks is completely similar. The 
conductance is given by 

G, = §[5(K.-V^o:,)+5(^Gt,-K=)] , (6) 

where Go is the conductance of the dot-reservoir barrier, 
Z = 2 + e""(^<^-^G:,) + e-"(^G..-^G) is ^j^g partition 
function with a — 2CgEc/ ek^T ^ and g{x) ~ Qa;/(e"^ — 

The calculated conductances for P and AP magnetiza- 
tion alignments as functions of normalized gate voltage 
for several temperatures are shown in Figs. 2a and 2b, 
respectively. As expected, the position of the resonant 




FIG. 2: Conductances for antiparallel magnetization align- 
ment (a) and parallel magnetization alignment (b), and TMR 
ratio as a function of normalized gate voltage eVo/Ec (c). 
Here Co = 0, Co = Cg, and kp = 2kap = Q.lEc- 

peaks depends on the magnetization alignment. The dis- 
tance between the split peaks is greater in the P align- 
ment than in the AP alignment, the resulting TMR ratio 
plotted in Fig. 2c. On either side of the two peaks for 
Gap, the TMR is positive, exhibiting the same sign as 
that in normal FM/NM/FM junctions. At low tempera- 
tures, e.g., k^T — 0.01i?c, the positive TMR may reach 
about 98%, corresponding to Gp/Gap — 100. It forms a 
striking contrast to the much smaller ratio Kp/kap — 2 
used in the calculation. Actually, for gate voltage Vg 
on either side of the two peaks for Gp, one can obtain 
Gp/Gap — exp[(Kp — «;Ap)/fcBT] , increasing exponentially 
with the temperature lowered. More interestingly, in the 
region between the two peaks of Gap, the TMR is not 
only very large in magnitude but also inverted in sign. 
For gate voltage in this region, the conductance ratio 
Gp/G AP — exp[— (Kp — Kap)/A:bT] decreases exponentially 
with the temperature lowered. With increasing temper- 
ature, both positive TMR and negative TMR decrease 
in magnitude. As energy k^T is increased so as to be 
comparable with the energy separation (kp — Kap) of the 
left (or right) two Coulomb peaks for Gap and Gp, the 
TMR ratio becomes very small. 

In Fig. 3, the conductance ratio Gp/Gap is plotted as 
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FIG. 3: Conductance ratio Gp /Gap as a function of nor- 
malized temperature kuT/Ec for different Kp and Cr>. Here 
Go — Gg and kp = 2k, ap- The gate voltage is Vg — V^ p 
for the upper three curves, and Vg = (VJ^p + Vq p)/2 for the 
lower three curves. 

a function of temperature for different interdot electro- 
static couplings Cd and interdot tunnel coupling . We 
note that Gp/Gap > 1 corresponds to positive TMR, and 
Gp/Gap < 1 corresponds to negative TMR. So far, we 
have confined ourselves to the simple case of zero inter- 
dot capacitor Gq. The interdot capacitor represents the 
effect of the excess charges in one dot on the electrical po- 
tential of the other dot. According to Eq. (O, a nonzero 
Co mainly enhances the splitting of the Coulomb peaks. 
The two dotted lines in Fig. 3 are obtained for nonzero 
Gd. By comparing them with the solid curves, we find 
that an increased interdot capacitance does not change 
the TMR significantly. 

Finally, we wish to make a discussion on experimental 
realization of the present theory. With today's fabrica- 
tion technology, it is possible to connect two FM dots by 
using short narrow conductor (quantum wire) [sJl, and 
keep the FM dots sufficiently separated in space. In this 
manner, relatively strong electron hopping coupling can 
be ensured, and in the meantime magnetic dipolar inter- 
action can be avoided between the FM dots. As a result, 
the RKKY interaction caused by electron hopping be- 
tween the dots may induce an AF alignment, if the dots 
are suitably distanced. The P alignment can be achieved 
simply by applying an external magnetic field. As a sim- 
plified model of the real systems, we have omitted the 
complexities such as spin-flip scattering and electron vir- 
tual tunneling in the Coulomb blockade regime. These 
effects could reduce the TMR measured experimentally. 
However, it is believed that they will not change the re- 
sults qualitatively. 

In summary, we have investigated the charge fluctua- 
tions and conductance peak splitting in magnetic double 
quantum dots. The characteristic energy of the Coulomb 
peak splitting is shown to depend on the magnetization 
alignment. Therefore, the change in the magnetization 
alignment causes very large TMR effect. The TMR ra- 
tio may be either positive or negative, depending on the 



value of gate voltage. 
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